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Quarkonium Production through Hard Comover Scattering. 
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We extend to large transverse momentum px an approach to quarkonium hadroproduction previ- 
ously suggested for low p±. The dynamics involves a perturbative rescattering of the heavy quark 
pair off a comoving color field which originates from gluon radiation prior to the heavy quark pro- 
duction vertex. Assuming simple properties of the comoving field we find the rescattering scenario 
to be in reasonable agreement with data. At large px, V' is predicted to be unpolarized, and xi 
production is favoured compared to X2- We predict the Xi polarization to be transverse at low px, 
and to get a longitudinal component at large px- 



I. INTRODUCTION 

Quarkonium production is a rich domain of QCD 
which involves several hardness scales, namely the trans- 
verse momentum p± of the heavy quark pair, the quark 
mass TO and the Bohr momentum of the heavy quarks 
in quarkonium ttgTO. Furthermore, bremsstrahlung glu- 
ons are produced along with the pair at an intermediate 
hardness scale. The proximity of the open flavour thresh- 
old makes the production of a quarkonium state highly 
sensitive to its production environment. Relatively soft 
rescatterings of the quark pair with surrounding fields 
can make or break the bound state. 

Quarkonium production may be contrasted with open 
heavy flavour production, for which the total cross sec- 
tion is unaffected by late rescattering. This insensitivity 
is technically expressed through the QCD factorization 
theorem 0] , which allows the cross section to be written 
in terms of process- independent parton distributions. An 
analogous factorization of hard and soft physics does not 
apply to quarkonium production rates, which constitute 
a small fraction of the total open flavour cross section. 
This makes quarkonium physics a challenging and poten- 
tially rewarding field which can teach us new aspects of 
the dynamics of hard processes. 

Attempts to ignore rescattering effects in quarkonium 
production have met with success in photoproduction 
but failed in hadroproduction. The PQCD prediction 



for quarkonium production without rescattering, at low- 
est order in the quark pair relative velocity v (the 'Color 
Singlet Model', CSM |,|), accounts well for the HERA 
data on 7p — > J/^X in the photon fragmentation region 
Q. On the contrary, the CSM underestimates the Teva- 
tron pp data on direct J/^ and ■0' production at large 
p± by a factor 50 Similar discrepancies are also 
observed in hadroproduction at low p± 0] . 

The CSM amplitudes contain infrared divergencies, 
which for P-wave states appear already at lowest or- 
der in as- The divergencies cancel systematically be- 
tween terms of different orders in w, as shown by the 
'NRQCD' 10] expansion of the amplitude, fn the absence 
of rescattering effects NRQCD is thus a theoretically vi- 
able framework for calculating quarkonium production in 
QCD. 

The NRQCD terms of higher order in v correspond 
to relativistic effects in the quarkonium wave function, 
and are related to its higher Fock components such as 
\QQg). The magnitude of the relativistic corrections is 
poorly known, and higher powers of v may be compen- 
sated by lower powers of as at the hard vertex. Hence one 
may consider the possibility that the observed discrepan- 
cies between the CSM and hadroproduction data are due 
to higher order contributions in the NRQCD expansion 
(the 'Color Octet Model', COM g). This impHes, how- 
ever, new contributions also in photoproduction which 
tends to lead to an overestimate of the data The 
phenomenological difficulties of a COM approach [PlJTc 
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have recently been made more acute by th e (p reliminary) 
CDF data on i/)' polarization at large p±_ fllf . The COM 
prediction of a transverse ip' polarization, considered as 
a crucial test of the COM approach, fails by more than 
3 standard deviations 

It would thus appear that rescattering is important 
in quarkonium hadroproduction. Since gluons, unlike 
photons, carry a comoving color field one expects more 
rescattering in hadroproduction than in photoproduc- 
tion. In Q (referred to as I in the following) a scenario 
involving a perturbative rescattering of QQ pairs pro- 
duced at low p± ^ m was studied. A simple form of the 
comoving field allowed several of the puzzling features of 
the data to be understood. 

Here we shall extend the rescattering picture of I to 
quarkonium production at high p±. We refer to I for a 
detailed discussion of the rescattering mechanism and its 
phenomenological justification. In section |l^ we present 
the low p± calculation in a more compact form than in I 
and recall the results which were obtained. In section [II 



we propose a generalization of the low pj_ approach to 
high p±, and evaluate S'-wave and P-wave quarkonium 
production amplitudes. All results obtained within the 
rescattering picture (at low and highpj_) are summarized 
in section IV, together with a discussion of our assump- 
tions. 
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FIG. 1. Scenario for the creation of a comoving color field 
in low p± hadroproduction. 

B. The model 



II. LOW P± QUARKONIUM PRODUCTION 
A. Physical picture 

The rescattering mechanism of quarkonium hadropro- 
duction proposed in I is based on an interaction of the 
heavy quark pair with a comoving color field. Let us re- 
call why a color field comoving with the QQ pair may be 
created in hadroproduction, but not in photoproduction. 

In the QED process e+e~ near threshold 

(Fig. la), the incoming electrons are surrounded by an 
electromagnetic (photon) field. As the electrons annihi- 
late, these fields pass through each other without inter- 
acting (except via ©(aem) processes). The pair is 
thus created in a field- free environment. 

In the analogous QCD process gg QQ (Fig. lb), the 
color fields associated with the color charge of the collid- 
ing gluons on the contrary interact strongly. Such multi- 
ple interactions may create a remnant color field with low 
rapidity components in the QQ rest frame. Rescattering 
may occur between this remnant field and the QQ pair. 

In photoproduction, jg QQ (Fig. Ic), the (unre- 
solved) photon carries no radiation field and the situation 
is analogous to QED: no remnant field is comoving with 
the QQ pair. Rescattering thus occurs only in hadropro- 
duction, and may explain some of the observed 'anoma- 
lies' of quarkonium production. It is absent in photopro- 
duction (except for resolved photon contributions), where 
the CSM predictions in fact are quite satisfactory. 



The effects on low p± quarkonium hadroproduction of 
a rescattering of the QQ pair off a comoving color field 
were studied in I in the framework of a simple model. 
Before extending this scenario to quarkonium production 
at high p± we summarize the main results of I. 

The model applies to high energy production, where 
the quarkonium carries a moderate fraction xp of the 
hadron beam momentum and has p± ^ to, where to is 
the heavy quark mass. The amplitude for quarkonium 
production in the bound state rest frame shown in Fig. 2, 

L~,Sz XX, aa 

(1) 

is a convolution of the gg — > QQ color octet wave function 
$[^1 , the kernel TZ describing the rescattering between the 
pair and the comoving color field, and the quarkonium 
wave function 'if. 

The gg QQ process occurs on a proper time scale of 
order m~^, whereas the time scale for quarkonium bound 
state formation is of order (a^TO)"^. The rescattering is 
assumed to be characterized by the 'semi-hard' scale fi 
of DGLAP gluon radiation 



ttsTO <C /i <C TO 



(2) 



As a consequence, the rescattering is well separated in 
time from both the QQ pair creation and the bound state 
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FIG. 2. The low quarkonium production amplitude in- 
duced by a rescattering on the comoving field F (two dia- 
grams). From Ref. j^. 



formation. The heavy quarks thus propagate nearly on- 
shell both before and after the rescattering. In Eq. (|l|) 
the hehcities and momenta of the quarks (in the quarko- 
nium rest frame) before the rescattering are denoted A, A 
and PtP' . After the rescattering of momentum transfer € 
the same quantities are denoted by a, a and q, —q. The 
calculation of I can be summarized by expressing every 
factor of (|l) as a contraction between the spinors xa and 



xL with x\ = (1 0), X- = (0 1), 

+ i27:)H'{p - q) 51 xU T^'^i^^ -1) X--a (3b) 



(3a) 



*<Tff (9) = -^J^X-s(i{Sz) Xo 



(3c) 



In TZ'^{i, q) and TZ'^{i, q), q stands for the momentum of 
the scattered quark or antiquark after the scattering. Eq. 
(^) holds for S — 1 quarkonia, e{Sz) being the bound 
state spin polarization vector defined in the rest frame as 



(4a) 
(4b) 



e(±l) = (Tl,-i,0)/V2 
e(0) = (0,0,1) 



We use the notation a = a - cr for a vector a, with Ufe the 
Pauli matrices. The amplitude (^ takes a form which 
can be readily inferred from Fig. 2, 

d'q nLM) 



M{S=l) = Y: {LLz;SSM) J A!^3 ^ 

L/z ,Sz 

-^T^, {Tt [eiSzrnQi£,q)m{q-£,-q) 



Tr 



n'^{£,-q)eiSzr¥'\q-q - £)]] (5) 



where Tr denotes the trace of 2 x 2 matrices. For con- 
venience we separate the rescattering color factor from 
the kernel TZ. The color indices of the rescattering gluon 
and of the quark and antiquark before the rescattering 
are denoted by 6, i, j. The general form ^ will also be 
used in the case of high pj_ quarkonium production via 
gluon fragmentation in the next section. 



In Eq. (g|) the zeroth and first order terms in q corre- 
spond to S and P-wave production, respectively. Recall- 
ing that \£\ ^ fi ra we work only to first order in the 
small quantities t/ni, q/m. The rescattering kernel then 
reads [01 



nQQ{£,q) = 



g 

2m 



TiK + G 



(6a) 



K{£, q) = 2mT°{i^,£) + T{i^,£) ■ {£ - 2q) (6b) 
G{£)=T{£°,£)x£ (6c) 

The color field comoving with the QQ pair is modelled as 
an external source T^{£), the form of which is not known 
a priori. The hardness scale /i of F is given by Eq. (^) . 

The wave function $[^1 of the quark pair produced in 
99 ^ QQ is § 



2m 



2m 



[ei ■ 62 Sp^as + ei • (5p (=2 + 62 • dp ii] 



(7) 



with Sp ~ p p', ai, a2 the color indices of the incoming 
gluons and — e(AA z = 1, 2 their polarization vectors. 
Using (|), (|) and (Q) we recover Eq. (18) of I for ^Si 
productionrn 



daia2b 2g^Ro 
4^/3 V2TTm3 



ei -£5^ 



62 -€(5^; 



} (8) 



Here we used for S'-wave bound states 



d^q 

(27r)3 



Rq 



(9) 



where Rq is the wave function at the origin. Summing 
over oi, a2, b and Ai, A2 yields 



,_iN^~A)iN^-l)4g^Rl 
3T:m^ 

(Sz = 0) 



167V, 

|G"|' + (r)2|F"(£)|' 
i|G^|Vi(£^)2|F0(£)| 



(10) 



-1) 



-Rc [iGyG"] 
where no assumption has so far been made on F^. 



*In Eq. (18) of I, the gluon propagator i/£'^ should be in- 
cluded in F^ and the color factor D'^ should be multiplied by 
51 5j . These details are not important for the conclusions of 
L 



3 



C. Results and predictions 

S-wave states are observed to be unpolarized at low 
p±_ p3|-p^. Comparing to Eq. (p^), this suggests that 
the gluons in F are dominantly transverse (|r°| <C |r|) 
and isotropically distributed. Hence we use the ansatz 



[{£) with A = ±1 



-Avw (11) 

where e^i{t} denotes a transverse polarization vector for 
the gluon £. The isotropy assumption enters in the inte- 
gral over £ after the sum over A is performed. In ( |l0|) we 
obtain 



A 



1^ 



lie 







(12) 



The treatment of P-wave production is similar. For low 
p± production the results can be summarized as follows: 

• Direct 5-wave production is unpolarized. 

• P-wave states are produced in the ratio 

CFdiriXl) / '^dir{X2) = 3/5. 

• Direct xi production is transversely polarized. 

• X2 is produced only with Jz = and — ±1 in 
the ratio adir{x2, Jz = 0)/(Tdir{x2, Jz = ±1) = 2/3. 

As noted in I, the decays of "^Pj states produced via the 
rescattering mechanism induce a longitudinal J/^ polar- 
ization. Since the J/tp is observed to be nearly unpolar- 
ized Jl3| , p^ , this indicates that also the CSM mechanism 
contributes to X2 production. CSM produced X2 states 
have Jz = ±2 only and decay to transversely polarized 
J/ip's. For acsMiX2) - crrescatt{X2) , an overall J/ ■(/; po- 
larization consistent with the data is obtained. Hence 
we also get a lower Xi/X2 ratio, adirixi) / o-dir{X2) 0.3, 
compatible with the data [l^ . 

It is interesting to observe that (low p±) P-wave pro- 
duction via rescattering is actually independent of F°. 
This is because only the S = 0, L = 1 part of the QQ 
wave function (|^) (the term proportional to faia2c) con- 
tributes in this case. Thus a spin- flip rescattering (trans- 
verse gluon exchange) is required to form a ^Pj state. 



the quark pair, and we assume that rescatterings occur 
between this field and the pair. The field comoving with 
a high p± QQ pair is a priori different from the one at 
low p± . 




FIG. 3. High p± quarkonium production via gluon frag- 
mentation. The QQ pair rescatters on the comoving color 
field Tf before forming quarkonium. 

A schematic view of our model is shown in Fig. 3 for 
'^Pj quarkonia (which require a single rescattering). A 
high p± transverse gluon (typically produced in a process 
like gg — > gg, which does not concern us here) fragments 
into a heavy QQ pair. The pair rescatters from an exter- 
nal field Tf oi hardness /i / similar to the scale fi of Eq. (|^) 
and forms a quarkonium bound state. The amplitude of 
the fragmentation process is a convolution, analogous to 



(|1|), of the quark pair wave function ^'^j^^j the reseat- 



Eq. 

tering kernel TZ)^x ,jg{£,p,q) given by Eqs. (3b) and 

and the quarkonium wave function -^^f'iq) of Eq. (g). 
As in the low p± situation, the heavy quarks are assumed 
to be on-shell before and after the rescattering. 

roi 

The wave function $ , , t of the pair in the quarkonimir 
rest frame is 



/,AA 



xl^T X^-x = 9Tt,uxipyiXg)v-^ip') (13a) 
^fip,p') = -2mgT5 {e(A,) + ^ e(A,) ^] (13b) 

Here a, i, j are the color indices of the gluon and quarks 
and e(Ag) is the polarization vector of the fragmenting 
gluon. For the range p± ^ m that we are considering the 
gluon is quasi-real and hence transversely polarized (Ag = 
±1) relative to its direction in the laboratory frame. An 
expansi on in the small ratios p/m, p' /m has been made 
in Eq. ( pb| ). 



B. Production of Si states 



III. QUARKONIUM PRODUCTION AT HIGH Px 

A. Scenario for rescattering 

Quarkonium production at high p± ^ m pro- 
ceeds dominantly through the fragmentation of quasi- 
transverse gluons |l^]. DGLAP radiation from the frag- 
menting gluon gives rise to a color field comoving with 



Two rescatterings are needed to produce states 
via gluon fragmentation, see Fig. 4. The fragmentation 
amplitude M / can thus be expressed as a convolution 
similar to (||) , with two rescattering factors TZ. Only the 
first term in the QQ wave function ( |13b| ), correspond- 
ing to a gluon fragmenting into a QQ pair with L — 0, 
actually contributes to S'-wave production. 
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FIG. 4. Highpx^5'i production in the rescattering scenario 
(four Feynman diagrams). 

We find for the '^S'l production amplitude 
3 c. o ^ -rngRo 



X { Tr [a 62 ^^i] Tr 



+ Tr[a62;^i]Tr [^'^(.ga, 0)e(5,)*7^'3(£l, 0), 

+ Tr[a6i62]Tr [e(5,)*^Q(.62, 0)e(A<,)^'3(£i, 0) 

+ Tr[a6i62]Tr [u"^ i£2,0)e{S,y e{Xg)n^ {£1, -£2)]} 

(14) 

where bi, 62 are the color indices of the two rescatter- 
ing gluons and Tr [a bi 62] is a shorthand notation for 
rp^. |-y ay 61^62 ]_ xhe four terms in (|l|) correspond to 
the four Feynman diagrams implied in Fig. 4. Using (|^) 
leads to the surprisingly simple result 



2V6^ 



Gi-e(Ag)G2-e(5,) 



(15) 

where Gi — T f{£i) x ii. Note that ( [l5| ) is independent of 
Fj. Only transverse gluons in the comoving gluon field 
contribute to ^Si quarkonium production at high p±. 

The expression (^) trivially satisfies the gauge in- 
variance requirement Mf(r'^{£i) — > if) = 0. Since 
the quarks are treated as being on-shell, each Feyn- 
man diagram contributing to the production amplitude 
separately satisfies current conservation. However, in 
order to keep track of gauge invariance when combin- 
ing the different kernels TZ, one must pay attention to 
their ^"-dependence. For instance the value of ii in 
n'^''^ {£1,-12) differs from that in TZ^''^ {£i,0) by an 
amount = -ii ■ £2/171. Now r'^{i'(,£i) is related to 
r^(£; + 5i°i,£i) through a boost of velocity (3 = -£2/^1, 
|/3| ^ 1. One realizes that the vector T fi^f' ,£) (and thus 
G) is insensitive to this boost (at the level of our ap- 
proximation) whereas the difference between T'^j-{t\,£i) 
5£i,£i) is such that we get 



and rO(£0 



nQ-^{£i,~£2)~nQ-^{£i,o) 

which gives the simple form (^5|) for A4 /. 



(16) 



In the S'-wave production amplitude ( p^ the depen- 
dence in Xg and Sz is factorized. No correlation appears 
between the fragmenting gluon and quarkonium polariza- 
tions. Consequently, assuming only that Fj is isotropi- 
cally distributed in the (comoving) QQ pair rest frame 
we predict that directly produced '^Si quarkonia are un- 
polarized at high p± . Contrary to the low p± situation, 
we do not now need to assume the dominance of trans- 
verse gluons in F/. The preliminary CDF data [|ll| on 
high p± ip' production seems to prefer a longitudinal po- 
larization, but the statistics is insufficient for a definite 
conclusion. 

The high p± cross section ratio a{ip')/a{J/ip) appears 
to be larger |^ than the nearly universal ratio measured 
in low p± hadroproduction ||l8|-^ , 



(Jdir{J/i^) 



0.4 > 



crdtr{J/ip) 



0.24 



P± 



(17) 



This is in qualitative agreement with our result that the 
rescattering gluons are transverse and thus couple to the 
spin of the individual quarks. In low p± production the 
target gluon is longitudinal (in the target rest frame) and 
its coupling is proportional to the QQ color dipole size 
r±. The convolution (|l]) then probes the bound state 
wave function in relatively larger configurations, where 
the node in the -0' wave function tends to decrease its con- 
tribution ||2l|,|2^ . This systematics carries over to diffrac- 
tive photoproduction, where two longitudinal gluon ex- 
changes give a factor , and the corresponding ratio (|lT 
is measured to be ~ 0.15 |2^. 



C. Production of ^Pj states 

In the case of high p± ^Pj production, only one scat- 
tering is needed (see Fig. 3). The production amplitude 
is given by the expression (^) with replaced by 
of Eq. (|l^) . We expand the curly bracket of (||) at q ^ 
and use for P-wave states 



{2n) 



= 



(27r)3 



(1'^*ilX(1) 



ATTm 



i?i'e(i. 



(18a) 
(18b) 



where Ri' is the derivative of the P-wave function at the 
origin. In the small q expansion the linear terms in q 
cannot arise from the rescattering kernel TZ because of 
the gauge invariance property mentioned in the preced- 
ing subsection. Recalhng ( [l6| ) with £2 replaced by —q, 
one indeed sees that 'RP{£,q) is actually independent of 
q. This means that the final orbital angular momentum 
of the QQ pair L = 1 must be fixed before the scattering, 
via the L — I part of the g — > QQ wave function l>j'. 
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This statement remains true for any numb er of rescatter- 
ings. Thus only the second term of (|l3b|) con tributes to 
P-wave production. Inserting (|6^) and (13b) in (||) we 
arrive at 

MfCPj, J,) = ^^^£i^JbM"^S"^''JeXAg)HJ (19a) 
A^"^ J2{LL,;SS,\JJ,)e''{L,re^{S,r (19b) 

B«/3.'J ^ SfS^i + SfS^ - d^Sf (19c) 
H = i^Tf + 2mT"f£ (19d) 

where B"^-'-^ arises from the trace of four Pauh matrices. 
Using the form of the tensor A"^ given in p3] we get 



Mf{^Pj,J.) 



e{\g) ■ H (J = 0) 

e(Ag) • (H X e(J,)*) (J = 1) 
(-^ = 2) 



(20) 



The amphtude vanishes for J = 2 since in this case A"^ is 
symmetric and A'^ — 0. ExperimentaUy (j{x2) — "'(xi) 
at large [ p5[ . The cross section for X2 production 
should be generated either by more than one rescattering 
in our scenario or by the CSM (or COM) mechanism of 
gluon emission. 

To conclude this section we give our prediction for the 
Xi polarization at large p± and discuss its dependence 
on the rescattering field T j. 

The high p_L fragmenting gluon being quasi transverse 
we have from (pO|) 



|Xf (J,)|^oc |(Hxe(J,)*)i|^ = 



|HxP (J. 



--0) 
±1) 



(21) 

where the direction of the high p± gluon in the labora- 
tory frame is chosen as the reference axis. For a general 
field Tf, the xi polarization depends on the relative im- 
portance of t he two terms contributing to the vector H 
given in ( |l9d| ). However, assuming that £ is isotropically 
distributed, we find the same result for the xi polariza- 
tion parameter 



aiJ, = +1) - a(J, = 0) 



(22) 



a{J, ^ +1) + a{J, = 0) 3 
in the two following extreme cases: 

(i) Tf contains only transverse gluons (the first term 
of H dominates) and is modelled by the ansatz ( |ll| ) 
used at low p± . 

(ii) Tf is purely longitudinal (the second term of H 
dominates) . 

Thus, rather independently of the form of Fj, we predict 
the xi yield to be longitudinally polarized if the rescat- 
tering scenario dominates xi production at high p±. 



IV. SUMMARY AND DISCUSSION 

We presented an extension of the approach of I to high 
p± quarkonium production. As in our previous work, we 
assumed the presence of rescattering between the heavy 
quark pair and a comoving color field, and found that 
such a scenario qualitatively agrees with the low and high 
p± data. We summarize here our main results and pre- 
dictions. 

At low , the rescattering mechanism applies only to 
hadroproduction, not to photoproduction (in the photon 
fragmentation region). In order to agree with the ob- 
served non-polarization of low p± hadroproduction, 
we assume the field F created at low p± to contain dom- 
inantly transversely polarized gluons. The cr(xi)/(T(x2) 
ratio is found to be compatible with the measured value 
and our main predictions at low pj_ are: 

(i) Direct xi production is transverse. 

(ii) (Tdirix2, Jz = 0)/(Tdjr(x2, Jz = ±1) = 2/3. In addi- 
tion, we expect a significant CSM (no rescattering) 
contribution to (t(x2), having — ±2. 

The comoving field F f arising from high p± gluon frag- 
mentation is a priori different from the low pj_ field F. 
Our scenario is identical in high p± hadroproduction and 
high p± photoproduction, the QQ pair originating from 
gluon fragmentation in both cases. Our predictions at 
large p± are: 

(iii) Direct '^S'l production is unpolarized. 

(iv) Direct xi production contains a longitudinal com- 
ponent (Eq. (H)). 

(v) X2 is not produced via (a single) rescattering. 

The prediction (iii) is consistent with the present CDF 
data on tp' polarization pT[ |, but higher statistics is 
needed for a definite conclusion. Only transverse gluon 
rescattering contributes to '^Si production. This is in 
qualitative agreement with the measured hadroproduc- 
tion ratio a{'ip')/adiriJ/i^) being smaller at low than at 
high p± (cf. Eq. (17)), and being even smaller in diffrac- 
tive photoproduction. 

All our results and predictions apply equally to the 
charmonium and bottomonium families. 

We used several simplifying assumptions which we now 
discuss. 

• We supposed the rescattering color field to be 
isotropically distributed in the QQ rest frame. Our 
predictions on ratios and relative polarization rates 
depend on this assumption. 

• We assumed the rescatterings to be perturbative. 
Thus only the minimal number of rescatterings re- 
quired to produce a given bound state was consid- 
ered. On the other hand the rescattering probabil- 
ity should be large enough for this mechanism to 
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dominate the CSM S'-wave and (low p±) xi rates. 
Whether such a compromise holds is not obvious 
and will be studied in a future work. 

• It is unlikely that the production amplitude is sen- 
sitive only to the quarkonium wave function (or its 
derivative) at the origin, at least for charmonium. 
This is particularly so for P-waves (see the discus- 
sion of (t(x2)/o'( J/'0) in I) and when 'dipole' factors 
of r± enhance larger wave function configurations 
0. 
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